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Abstract: 

Iterative  solutions  of  coupled  parabolic  wave  equations  are  examined  to 
determine  the  validity  of  applying  this  method  to  predict  low  frequency  ocean 
acoustic  reverberation.  Solutions  for  backscattering  from  a  random  rough 
interface  between  two  media  of  differing  wave  number  are  obtained  in  the  form 
of  an  expansion  in  powers  of  the  rough  interface  excursion  from  its  mean 
value.  This  expansion  is  found  to  differ  significantly  from  the  corresponding 
expansion  obtained  by  applying  perturbation  theory  to  the  full  elliptic  wave 
equation. 
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INTRODUCTION 


The  parabolic  approximation  has  been  successfully  applied  to 
electromagnetic  propagation1,  non-linear  acoustic  phenomena2,  and  underwater 
acoustic  transmission3.  Claerbout*  appears  to  first  have  considered  iterative 
solution  of  a  coupled  pair  of  parabolic  equations  to  obtain  the  backscattered 
field  applicable  to  seismic  reflection  profiling.  Because  the  parabolic 
approximation  provides  a  computationally  efficient  numerical  method  for 
solving  wave  propagation  and  scattering  problems,  its  applicability  to 
underwater  acoustic  backscatter  is  well  worth  investigating. 

In  this  article,  a  two-dimensional  pair  of  coupled  parabolic  equations  is 
considered.  The  solution  to  this  pair  of  equations  is  obtained  for 
backscattering  from  a  random  rough  interface  between  two  media  of  the  same 
density  but  differing  indices  of  refraction.  To  facilitate  comparison  with 
other  scattering  theories,  results  are  obtained  for  a  plane  wave  incident  on  a 
rough  surface  area  in  the  far  field  of  a  receiver,  and  ensemble  averaged. 

Perturbation  theory,  which  agrees  well  with  exact  Monte-Carlo  simulations 
of  scattering  from  pressure  release3  and  fluid-solid6  interfaces  of  moderate 
roughness,  is  used  for  comparison.  For  the  lowest  order  scattering  which  is 
due  to  Bragg  diffraction,  the  parabolic  and  perturbation  results  are  found  to 
differ  significantly.  Higher  order  scattering  terms  are  also  evaluated  with 
similar  findings. 

Section  I  of  this  article  presents  the  general  formulation  of  the  problem 
in  which  the  elliptic  wave  equation  is  split  to  obtain  a  coupled  pair  of 
parabolic  equations.  In  Section  II,  this  coupled  pair  of  equations  is  solved 
to  obtain  the  backscattered  field.  Numerical  comparisons  between  perturbation 
theory  and  the  results  of  Section  II  are  presented  in  Section  III.  Section  IV 
summarizes  and  discusses  the  results. 


I .  GENERAL  FORMULATION 

The  two-dimensional  wave  equation  in  a  medium  having  a  variable  index  of 
refraction  k(x,z)  is  considered 


~r~r  +  +  k2  (x,  z)  v  =  0  , 

ox 2  dz2 


where  v  is  the  acoustic  pressure.  For  acoustic  energy  propagating  close  to 
the  x  direction,  it  is  convenient  to  rewrite  this  equation  in  matrix  form  as 

d  0  1]  C2-' 

Si  *  x  L.  0  ¥  ' 

dx)  V  7  dx) 

where  A2  -  k2(x,z)  +  d2/dz2.  Equation  (2)  may  be  split  by  a  method  due  to 
Corones7  to  obtain  a  pair  of  coupled  equations  for  the  transmitted  and 
reflected  fields,  v+  and  v.  respectively,  where 


■ T  m  ■ 


and  T  is  the  splitting  matrix.  Claerbout's4  result  follows  from  the  choice 


r  -  f 1 

I  1  *i/A ) 


Substitution  of  Eq.  (3)  into  Eq.  (2),  then  yields  the  following  pair  of 
equations 


^  (v.-vJ-iAv. 


3 


Equations  (5)  and  (6)  are  exact;  no  approximations  have  been  made  up  to 
this  point.  However,  to  obtain  an  efficient  numerical  solution  of  this  pair 
of  equations,  it  is  usually  assumed  that  the  reflected  field  is  small  and  may 
be  neglected  in  the  solution  of  Eq.  (5),  or  that 


dv, 

dx 


dA 

dx 


v,  +  iAv ;  , 


(7) 


which  with  a  suitable  approximation  for  A,  may  be  solved  for  the  transmitted 
field.  With  the  transmitted  field  known,  the  reflected  field  may  be  found  by 
solving  Eq.  (6). 

This  procedure  for  numerically  solving  Eqs.  (6)  and  (7)  is  illustrated 
schematically  in  Figure  1.  A  parabolic  equation  is  numerically  solved  for  the 
transmitted  field  v+  out  to  some  maximum  range  of  interest.  Values  of  v+  are 
stored  at  each  point  where  3A/3x  is  nonvanishing.  Equation  (6)  is  then  solved 
for  the  reflected  field  with  the  initial  condition  that  v.  -  0  at  the  maximum 
range  of  interest.  If  desired,  with  v.  known,  one  may  return  to  Eq.  (5)  to 
refine  the  estimates  of  v+,  and  iterate  again  to  better  estimate  v. .  This 
method  is  certainly  numerically  tractable  for  scattering  from  rough 
interfaces,  since  only  the  transmitted  and  reflected  fields  on  the  interfaces 
need  be  stored. 


II.  BACKSCATTER  FROM  ROUGH  INTERFACES 
A.  The  Transmitted  Field 

To  determine  if  the  method  described  in  Section  I  does  indeed  provide  a 
good  estimate  of  the  reverberant  field,  we  consider  scattering  at  an  interface 
between  two  media  having  the  same  density  but  differing  indices  of  refraction. 
For  simplicity,  an  incident  plane  wave  is  assumed  as  illustrated  in  Figure  2. 
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To  obtain  an  expression  for  the  transmitted  field,  the  operator  A  is 
approximated  as 


A  =  k(x,  z)  + 


i  a2 

2k0  dz2  ' 


where  k0  is  a  reference  wave  number.  In  Appendix  A,  it  is  shown  that  the 
first  term  on  the  right-hand  side  of  Eq.  (7'  contributes  negligibly  to 
backscattering.  With  the  neglect  of  this  term,  the  transmitted  field  above 

the  interface  v*15  then  obeys 


(d_iK  i  JL]  v<»  .  o 

^dx  2k0  dz2) 


and  below  the  interface  v,  obeys 


( J. .  iK  -  JL  JL) 

\dx  2 k0  dz2) 


v<21  =  0  , 


where  Kj  and  Kz  are  respectively,  the  (constant)  indices  of  refraction  above 
and  below  the  interface. 

For  an  incident  plane  wave 


vin  -  exp  [i(kx-kjZ)]  , 


v*1J  and  vf2>  are  given  by 


»  exp  (ikx)  [exp {-ikxz)  +  Rex.p  (iklz)] 

m 

♦  J  dk'exp  [i  (k'x  *  k[z)  ]  u^(k')  , 
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(13) 


vj2)  *  T  expli  (kx-k2z)  ] 

m 

+  j*  dk'exp  [i  (k'x  -  k'2z)  ]  u2(k')  , 
o 


where  kj  =  [2k0  (JCj  -  k)  ] 1/2  ,  and 
R  -  (k^kj/ik^kj  , 

T  =  2k1/(k1  +  k2)  . 

The  quantities  ux(k')  and  u2(k')  may  be  expanded  in  a  functional  series  in  the 
rough  surface  extension  f(x)  from  its  mean  value.  Following  Kuperman8,  one 
obtains 


u^k')  =  (ki-ki)  T  (ih(k'-k) 

ttt 

+  f  dk"  h(k'  -  k")  h(k"  -  k) 


k"-k"+-<***k*) 


(14) 


u,(k') 


(ki  -  k[)  T 


*m  I 

ih(k'-k)  *  f  dk"h {k'  -  k")  h(k" -k)  k'{  - kj' *  -(^  ~  *x) 


+  ,  , 


(15) 


where  h(k)  is  defined  by 


CU)  =  j  dkh(k)  exp(ikx)  . 


(16) 


B.  The  Backscattered  Field 

To  obtain  the  backscattered  field,  A  is  approximated  as  in  Eq.  (8)  and 
the  term  v_3A/3x  appearing  on  the  right-hand  side  of  Eq.  (6)  is  neglected  (see 
Appendix  B) ,  whence  v.(1)  obeys 
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where  the  Green's  function  G-11  is  given  by 

Gll)  =  -  [k0/2it  (x'  -  x)  ] 1/2  exp{-i  [jt/4  +  Kx  (x  -  x')  *  k0  (z  -  z')  2/2  (x  -  x')  ] }  ,  (^) 

for  (x  -  x')  <0  and  vanishes  otherwise.  The  geometry  of  Figure  2  will  be 
considered,  with  the  receiver  at  x  -  0,  z  -  z0,  and  the  scattering  region 
located  at  x'  -  L  ±  Ax/2,  where  L  »  Ax  and  Ax  »  X ,  the  acoustic  wave  length. 

In  this  case  gJ11  may  be  simplified 

G.11  -  -  (ic0/2icL)1/2exp{i  [<J>  *  Kxx' -  k0zQz' / L])  ,  (20> 

where  4>  -  ~x/U  +  k0z02/2L  and  k0z'2/2L  is  neglected,  since,  as  is  apparent 
below,  the  only  values  of  z'  that  contribute  are  on  the  rough  interface 
z'  -  C(x') . 

For  an  interface  between  two  media, 

k{x,z)  =  KxH[z  -  C(x)]  +  K2H[Z(x)  -  z)  ,  (21) 

where  H(x)  -  {  q,’  So  . 

With  the  definition  of  A,  and  Eq.  (21) 

=  i  A"2  k{x'z)  •  (22) 

The  scattered  field  may  then  be  found  by  substituting  Eqs .  (12),  (20)  and 
(22)  into  Eq.  (19) 


7 


viu  (0,  z0)  -  (K2-Kx)  f  dx'M*  4(J^>3  iSisl  e‘1>A-*v.[x'.C(*/)]  ,  <23> 

t-ix/2  ^ 


where  the  integration  over  the  rough  surface  has  been  replaced  with  an 
integration  along  the  x'  axis. 

By  expanding  the  exponentials  that  appear  in  v+  and  G.  as  a  power  series 
in  the  surface  excursion,  the  corresponding  series  for  v.  may  be  obtained.  In 
the  following,  we  can  approximate  k[x',f(x')]  A"2  as  1/KX  since  only  low 
grazing  angles  are  of  interest.  With  these  approximations 


L*  Ajf/2  f 

(0 ,z0)  =0  /  cix'  ifj£iex p(ijqx')  (a3 +  a2  +  a3  •*■  . . .)  , 

L-lx/2  °x 


where 


a3  =  (1  +  R)  exp(ikx')  , 

a3  *  (x1)  exp  (ilex')  [kxll-R)  +  ( k0z0/L )  (1+/?)] 

m 

+  iT  f dk'exp  ( ik'x1)  (k'2  -  k[)  h(k'  -  k)  , 


a3  =  -Q3lx/)  exp(ikx')  [  (k^Zo/L)  (1  -  R)  +  (k2  +  klzl/L 2)  (1  +J?)/2] 

+  Tf  dk'exp  (ikV)  (k2  -  k()  if  dk"  hlk'  -  k")  h  Ik"  -  k)  k('  -  k"  +  (k2  +  k2} 
-  ((*')  (k[  -kx)  h(k'-k)}  , 

-IK3-kx)  (  k0  V/2  ** 

and  ff  ■  — e»  . 


C.  Backscatterine  Strength 


It  is  evident  from  Eq.  (24)  that  v_<11  has  the  form 
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(26) 


v.11  *  wx  +  w2  *  w3  +  .  .  .  , 
where  the  subscript  on  w  denotes  its  order  in  the  rough  surface  excursion  from 
its  mean  value.  The  expectation  <  >  of  the  scattered  intensity  is  then  given 
by 

<vi1,v:w,>  =  <w1w1’>  +  <w2wj>  +  <w1w3*>  +  <WiWj>  +  .  .  .  .  (27) 

The  first  term  on  the  right-hand  side  of  this  expression  is  the  Bragg 
scattering  term  and  is  the  major  contributor  to  backscatter  for  slightly  rough 
surfaces.  The  remaining  terms  are  higher  order  corrections  which  become 
important  for  moderately  rough  surfaces.  For  very  rough  surfaces,  one  expects 
the  series  to  diverge  as  it  does  in  perturbation  theory. 

Substituting  Eq.  (16)  into  Eqs .  (24)  and  (25),  yields 

< Wj w j >  =  (A x/L)  [  (K2  -K2)  {Kr  +  k)  ]2k0  TT'W(K2  *  k)  /  [4  (2)1/zxf]  (28) 

where  f (x)  has  been  assumed  Gaussian  with 

<h(k)h'(k')>  =  w(k)  6(k-k')  , 
and  the  rms  surface  excursion  a  obeys 

+«■ 

o2  =  J  dkW(k)  . 

— m 

The  factor  of  (2)1/2  appearing  in  the  denominator  of  Eq.  (28)  arises  from 
considering  a  Gaussian  illumination  pattern  as  in  Reference  9. 

Before  continuing,  note  that  some  simplifications  result  if  monostatic 
backscatter  is  considered  and  k0  is  selected  equal  to  k  so  that  k0z0/L  -  kx. 

In  this  case 
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(29) 


<wJVj>  «  %nbxk\QQ*  J  dk'w(k')  W{2k  +  k') 

x  {4k2  -2k2[f(k')  +  f(k')]  ♦  [ic2  *  (k  +  k')2]  f(k')  f  Ik')}/ (2) 1/2  , 
where  f(k')  -  (k2*  -  kj/)/^  +  kx)  ,  and 

<w3Wi>  =  &nbxk2QQ'T’W(2k){-2klo3 

m 

+  rf  dk'w(k  +  k')  [ka(ka  -kt>  +  (£'-£')  (£[  +  k2  -  2kx)  ]}/  (2) 1/2  , 

-m 

where  Rj  *  [2k0  (Kj  -  k1  -  2k)  ] 1/2  . 

To  obtain  Eqs.  (29)  and  (30),  we  have  used 

<h(k)  h(k')  h(k")  h(k'")>  =  W(k)  W(k")  &(k  +  k')  &(k"  +  k'") 

*  W{k)  w(k')  [6  (k  ♦  k")  5  <k'  *  k'")  +  5  (k  ♦  k'")  5  (k'  *  k")  ]  , 


(30) 


and  h(k)  =  h"(-k)  .  It  has  also  been  assumed  that  Ax»Lc  the  rough  surface 
correlation  length,  and  that  Kx  *  k . 

Finally,  the  scattering  strength  SPE  may  be  obtained  from  <v.(1,vi1>’> 

SPB  =  (L/Ax)  <vlv  v^v’>  .  (31) 

It  is  worthwhile  remarking  that  the  dominant  contributions  to  Eqs.  (29)  and 
(30)  arise  from  the  first  terms  in  curly  brackets  on  the  right-hand  sides  of 
these  expressions.  The  remaining  terms  are  negligible  for  the  cases 
considered  in  the  following  section.  Physically  these  terms  and  the 
contributions  discussed  in  Appendix  A  correspond  to  multiple  forward  scatter 
and  shadowing . 
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III.  COMPARISONS  WITH  PERTURBATION  THEORY 
A.  Bragg  Scattering 

Analogous  to  the  development  of  Section  II  A  for  the  parabolic  equation 
(PE) ,  in  perturbation  theory  (PT) ,  the  pressure  p  in  the  upper  and  lower  media 
due  to  an  incident  plane  wave 

pln  =  exp  [iikx-i^z)]  ,  02) 

may  be  expressed  as 


and 


Pj  =  exp(ikx)  [exp(-iKiZ)  +  R0exp  (i^z] 
+  j  dk'exp  [i  {k'x  +  x[z)  ]  gx  {k')  , 


(33) 


p2  =  T0exp  [i  {kx  -  KjZ)  ] 

+  J  dk'exp  [i  {k'x  -  k£z)  ]  q2  {k')  , 


(34) 


where  xi  =  (Kli-k2)1,i  ,  and  R0  and  T0  are  the  Rayleigh  reflection  and 
transmission  coefficients. 

Again,  following  Kuperraan,  a  functional  series  for  the  qj  may  be  derived. 
The  backscattering  strength  is  then  found  from 

S„  =  K1sin0<g1  {-k)  g{  {-k)  >  . 

The  lowest  order  contribution  to  qx  is 

g1{-k)  =  -2i  k2  Ra  h(-2k)  , 

and  the  corresponding  backscattering  strength  is 

SpT (1,1)  =  4k|r0R0*  sin0  W{2k)  . 

The  analogous  Bragg  term  for  the  PE  case  can  be  written 


(35) 


(36) 


(37) 
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W(k  +  Kx)  . 


(38) 


SPg{  1,1) 


Ak3  ( ~  )  2  V*i  *  k)  3 

|ic2  +  iCi|2  4(2 ) 1/2 


On  comparing  Eqs .  (37)  and  (38),  it  is  evident  that  the  functional  dependence 
of  the  two  expressions  differs  at  low  grazing  angles  with  Eq.  (37)  varying 
with  the  fourth  power  of  the  grazing  angle  and  Eq.  (38)  with  the  second  power. 

To  illustrate  the  differences  in  the  PE  and  PT  Bragg  results,  we  consider 
a  surface  wave  number  spectrum  W(k)  obeying 


W(k)  =  c/k 3  ,  KL<\k\<Km  ,  (39> 

which  yields  a  result  that  is  independent  of  frequency  at  lowest  order. 

Figure  3  compares  the  PE  and  PT  predictions  for  c  —  0.001  and  K2/Ki  —  0.985 
corresponding  to  a  critical  angle  of  ten  degrees .  Although  the  two  curves  are 
similar,  differences  of  approximately  10  dB  occur  at  the  critical  angle  and  at 
very  low  grazing  angles.  Figures  4  and  5  compare  the  results  for  a  critical 
angle  of  five  degrees,  and  an  angle  of  intromission  of  five  degrees, 

K2/Ki  ”  0.996  and  1.004,  respectively.  Even  for  these  cases,  where  the 
scattering  is  weak,  significant  differences  between  the  two  predictions  are 
apparent . 

In  Figures  3,  4,  and  5,  k0  -  K:.  From  the  definition  of  kx  and  k2,  and 
Eq.  (38),  it  is  clear  that  k0  is  a  multiplicative  factor,  and  that  adjusting 
this  factor  within  reasonable  limits  will  not  significantly  alter  the  results. 
Thus,  the  disparity  between  the  PE  and  PT  Bragg  scattering  predictions  does 
not  appear  to  be  a  consequence  of  the  angular  propagation  error  inherent  in 
the  parabolic  approximation. 

It  is  evident  from  the  analysis  of  Section  II  that  each  operation  by  the 
operator  3A/3x  introduces  a  factor  of  f  in  the  solution  of  Eqs.  (5)  and  (6) 
for  backscattering  from  rough  interfaces.  Hence,  although  a  number  of 
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approximations  have  been  made  in  obtaining  Eq.  (38),  none  of  these 
approximations  affects  this  lowest  order,  or  Bragg  scattering  result.  It  is 
also  clear  that  repeated  iteration  of  Eqs.  (5)  and  (6),  as  described  in 
Section  I,  can  only  lead  to  terms  of  higher  order  in  the  surface  excursion. 

It  thus  appears  that  solution  of  the  full  coupled  pair  of  equation  is  needed 
to  accurately  model  Bragg  diffraction. 

B.  Higher  Order  Scattering 

For  moderately  rough  surfaces,  higher  order  terms  in  PT  lead  to  enhanced 
scattering.  The  second  order  contribution  to  qx  is 

g1(-ic)  *  f  dk'h(-k-k')  h{k'-k)  (x2  +kJ ~x!2)  ,  (4°) 

and  the  corresponding  backscattering  strength  is 

S„( 2,2)  =  lexfotfosine  jdk'w(k*k')  W(k-k')  |k( -k*  +  k2|2  .  (41) 

o 

The  comparable  PE  prediction  is 

SPB( 2,2)  =  (L/ Ax)  <Wjwi>  , 
where  <w2w2*>  is  given  by  Eq.  (29). 

To  compare  the  predictions  of  Eqs.  (41)  and  (42),  we  select  KL  -  0.01  Kx 
and  Km  -  2KX  in  Eq.  (39).  This  choice  provides  predictions  independent  of 
frequency  and  corresponds  to  an  rms  roughness  a  -  14  m  for  a  frequency  of 
200  Hz.  Figure  6  compares  the  (2,2)  PE  and  PT  results  for  a  ten  degree 
critical  angle.  As  in  Figure  3,  the  PE  predictions  exceed  the  PT  result  at 
very  low  grazing  angles  and  are  significantly  lower  near  the  critical  angle. 
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It  Is  also  of  Interest  to  examine  the  remaining  terms  that  contribute  to 
the  fourth  order  scattering  strength  prediction.  The  third  order  contribution 
of  qr  takes  the  form 


( -k)  *2  iK1j?0  f  dk'dk"  h(-k  - k ')  h  ( k'  -k")  h(k"~  k) 

(k'-kI)  -k>UL±M  +  kJJs!LJzL 

2  2  (Kj  Kj)  2  2 

_  +  (y HI  -  k"(k'  -  k«)  -  -  Kz( 


The  scattering  strength  is  then  given  by 


(43) 


-Kfc  (Ka  -k^) 

+  ( Kj  ♦  ) 


Spj.U.3)  -  4k?J?0J?o*  f/(-2Jc)  |  dk'W(k-k') 

—m 

-  (Kj-kJ)  (Rj-Ri)  ♦  (y*  -K^Kj  ♦  Kj'3)  ♦2K1(Kj-Ki)]  , 
where  R^  *  [K§  -  (k'  -  2k) 2] 1/3  .  The  corresponding  result  for  the  PE  case  is 


(44) 


SPB (1,3)  =  (L/Ax)  <w1w,’>  ,  (45) 

where  <w1w3*>  is  given  by  Eq.  (30). 

Figure  7  compares  the  full  four'_n  order  scattering  strength  predictions 
for  a  critical  angle  of  ten  degrees  for  c  -  0.001  and  the  same  choice  of  KL 
and  as  in  Figure  6.  For  PT,  the  sum  of  Eqs.  (37),  (41),  and  (44)  plus  its 
complex  conjugate  is  shown.  The  PE  result  is  the  sura  of  Eqs.  (38),  (42),  and 
(45)  plus  its  complex  conjugate.  In  the  PT  case,  the  (1,3)  plus  (3,1) 
contributions  are  positive  below  the  critical  angle,  and  negative  above  it. 
The  corresponding  PE  contributions  are  negative  at  all  grazing  angles,  hence 
the  increased  disparity  slightly  below  the  critical  angle. 
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IV.  DISCUSSION 

A  coupled  pair  of  parabolic  equations  has  been  solved  for  the  intensity 
backscattered  from  a  random  rough  interface.  The  result  has  been  expressed  in 
the  form  of  a  series  in  the  rough  interface  excursion  for  comparison  with  the 
results  of  perturbation  theory.  Even  the  lowest  order,  or  Bragg,  term  in  the 
parabolic  approximation  differs  significantly  from  the  perturbation  theory 
result.  Because  all  terms  of  the  first  order  have  been  retained  in  the 
iterative  solution  of  the  parabolic  equations,  this  disparity  can  only  be 
attributed  to  the  original  splitting  of  the  elliptic  wave  equation  and  neglect 
of  the  backscattered  field  in  the  initial  solution  for  the  transmitted  field. 
Higher  order  corrections  to  the  scattering  strength  follow  the  same  trends  as 
the  Bragg  terms,  with  significant  differences  between  the  parabolic  and 
perturbation  results. 

The  parabolic  approximation  offers  several  advantages  for  studying 
scattering  from  rough  interfaces:  it  is  computationally  efficient,  includes  a 
full  diffractive  treatment  of  shadowing  effects,  and  can  handle  focusing 
effects  at  the  scattering  surface.  For  these  reasons  it  is  useful  in  the 
study  of  forward  scatter  for  which  its  predictions  are  in  good  agreement  with 
perturbation  theory.  However,  due  to  the  disparities  noted  above,  iterative 
solutions  are  of  doubtful  utility  for  the  study  of  reverberation. 
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APPENDIX  A 


To  estimate  the  results  of  including  the  term  -A*13A/3x  v+/2  on  the 
right-hand  side  of  Eq.  (7)  consider 

v.  *  v;  *  . 

where  to  lowest  order,  from  Eq.  (12) 

vl  ■  Texp  (ikx)  , 

and  v.  is  treated  as  a  small  perturbation,  obeying 

SlgL  .  (A-D 

The  appropriate  Green's  function  for  the  solution  of  Eq.  (A-l)  is 

G.  =  [k(x,  z)  /2n  (x  -  x')  ] 1/2  exp{i[n/4  ♦  k(x,  z)  ( x-x ')  *  k0  (z  -  z')  2/2  [x  -  x')  ] }  , 

(A-2) 

for  0<x'<x.  To  proceed,  we  select  k0  -  Kx  and  approximate  k(x,z)  -  Kx. 

Solving  Eq.  (A-l)  via  Eq.  (A-2),  then  yields 

P.(x,  o  *  Q'T  [  dx' ■%-  (x-x')‘1/2  exp/iJqU'"-  (C  -  (')  2/2  (x  -  x')  ] }  .  <A'3> 

{  dx' 

where  Q'  -  (L)1/2Q.  Neglecting  the  second  factor  in  the  exponential  of 

Eq.  (A-3)  is  valid  for  small  surface  slopes.  Then  by  using  Eq.  (16)  for  f(x) 

and  making  the  change  of  variables  x'  -  x-q,  Eq.  (A-3)  takes  the  form 

X 

v,{x,0  *  Q'T  f  ik'dk' h(k')  explixiKj^*  k')]  j  exp(-igA')  .  (A-4) 

Extending  the  upper  limit  of  the  q  integral  to  ®,  then  yields 
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C)  *  Q'T  j  ik'dk'h(k')  exp[ix(Kx*  k')  *  iSn/ 4]  (it/2  |/c'|) 1/2  .  (A-5) 

— M 

where  S  -  1  for  k'<  0  and  S  -  -1  for  k'  >  0. 

To  lowest  order,  the  resultant  backscattered  field  v.(0 ,zQ)  is  given, 
from  Section  II  B,  by 

L*  Ax/2 

VA0.Z0)  =Q  f  dx'  iSi4iexp(i/r.x/)^[x/,((x/)J  .  (A-6) 

l-L,  2 

Substituting  Eq.  (A-5)  into  Eq.  (A-6),  yields 

£♦  Ax/ 2 

=  QQ'T  /  dx'  exp  (2ik1x/)  f  ik'dk'h(k')  exp  (ix/e'  ♦  iSn/4)  (rt/2  | Ar'l ) 1/2  . 

L-ix/2  °X  -m 

(A-7) 

On  comparing  Eq.  (A-7)  with  the  a2  term  of  Eq.  (24),  it  is  apparent  that  an 
additional  factor  of  Q'  appears  in  v. .  Because  |Kj  -  K2\/K1  «1,  in  the 
examples  treated,  the  contributions  arising  from  the  (7.  term  are  negligible. 
One  may  verify  this  computationally  by  evaluating  such  terms  as 

m 

<  =  \QQ'T\2n*bx I  dk'fUKi  -  k')  W(KX  *  k') 

0 

x  ((Xi  -k')2\Kx  ♦  *'!  +  (tfj+kVlAi  -k'\  +  2  (A?  -k>2)1'2H(K1  -k,)}/{ 2) 1/2  . 
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APPENDIX  B 


Following  the  method  of  Appendix  A,  the  result  of  retaining  the  term 
-A'13A/3x  v./2  on  the  right-hand  side  of  Eq.  (6)  may  be  estimated  by 
considering 

v.  -  vl*  V.  .  (B-l) 

where 

IT  =  j  dx'dz'A'1  v.(1>  (x',  z')  G.  (x  -  x! ,  z  -  z')  /2  ,  (B-2) 

and 

=  -j  dx'dz'A'1  ~  vl(x',  z')  G.  (x  -  x/,  z  -  z')  /2  .  (B-3) 

For  x'>x,  to  lowest  order 


:<x,C)  -  (?Tf 


—  d*'  ^rr-eicpif1(2x'-x) 
( x'-x)1' 2  dx'  1 


(B-4) 


where  kQ  «  and  k(x,z)  •  Kx .  Equation  (B-4)  may  be  integrated,  as  in 
Appendix  A,  yielding 


v_(x,C)  *  i  (x/2)  1/2Q't  J 


k'dk'h  (k')  [i(k'x*  K  x  +  sn/ 4)1  , 

\k'  +  2Kl\in 


(B-5) 


where  S  -  -1  for  (k'  +  2KX)  <  0,  and  S  -  +1  for  (k'  +  2K:)  >  0.  Substituting 
Eq.  (B-5)  into  Eq.  (B-3),  then  yields 


P.(0  ,z0) 


L’lx/2  / 

-i  (ji/2)  1/2qq't  f  dx'  - 

L-L/2  d* 


7  k'dk'h(k') 

im  \k'-2Kx\1/3 


x  exp  [ix'(2Kx  *  k')  *  iSn/A]  . 


(B-6) 


18 


As  an  example,  we  consider  the  expectation  <v_v!>  to  which  only  the  q2 


terra  of  Eq.  (24)  contributes 


<P_v!>  -  QQ*Q,T4bxni/2k1K1  j  dk'w(k')  W(2K1*kl) 

-m 

x  Sk'llK^  *  k'\x,i  exp(iSx/4)  , 

which  is  smaller  by  a  factor  of  approximately  (1  -  Kj/l^)  than  the 
contribution  of  Eq.  (29). 


(B-7) 
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Figure  6.  Higher  Order  (2,2)  Contributions  to  Backscatter  for  a  Critical  Angle  of  Ten  Degrees 
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Figure  7.  Fourth  Order  Backscattering  Strength  Predictions  for  a  Critical  Angle  of  Ten  Degrees. 


Distribution  List  for  Unclassified  ARL  Penn  State  TM  91-104,  entitled, 
"Backscattering  from  Rough  Interfaces  and  the  Parabolic  Approximation, 
by  S.  T.  McDaniel,  dated  12  April  1991. 


Scientific  Officer  Code  1125AO 
Office  of  Naval  Research 
800  North  Quincy  Street 
Arlington,  VA  22217-5000 

Attn:  M.  H.  Orr  -  Code  1125AO 
Copies  No.  1-3 


Administrative  Grants  Officer 
Office  of  Naval  Research 
Administrative  Contracting  Officer 
The  Ohio  State  Univ.  Research  Ctr. 
1314  Kinnear  Road 
Columbus,  OH  43212-1194 

Attn:  Resident  Representative  N66005 
Copy  No .  4 


Director 

Naval  Research  Laboratory 
Washington,  DC  20375 

Attn:  Code  2627 

Copies  No.  5-10 


Defense  Technical  Information  Center 
Building  5,  Cameron  Station 
Alexandria,  VA  22314 

Copies  No.  11-22 


Applied  Physics  Laboratory 
University  of  Washington 
1013  N.  E.  40th.  Street 
Seattle,  WA  98105 

Attn:  Eric  H.  Thorsos 
Copy  No.  25 


Office  of  Naval  Research 
800  North  Quincy  Street 
Arlington,  VA  22217 

Attn:  Ralph  Baer  -  Code  1125AO 
Copy  No.  26 

The  Pennsylvania  State  University 
Mechanical  Engineering  Department 
University  Park,  PA  16801 

Attn:  Alan  Pierce 
Copy  No.  27 


Applied  Research  Laboratory 
The  Pennsylvania  State  University 
Post  Office  Box  30 
State  College,  PA  16804 

Attn:  ARL  Library 
Copy  No.  28 


Naval  Underwater  Systems  Center 
New  London,  CT  06320 

Attn:  Ding  Lee  -  Code  3332 

Copy  No.  23 


Naval  Ocean  and  Atmospheric 
Research  Laboratory 
Stennis  Space  Center,  MS  39529-5004 


Attn:  Stanley  Chin  Bing  -  Code  221 
Copy  No.  24 


